ЗА 


УДК 004.93 


М.А. Алексеев, И.М. Удовик, К.А. Ахметшин 
ГВУЗ «Национальный горный университет», г. Днепропетровск, Украина 
Украина, 49005, г. Днепропетровск, пр. К. Маркса, 19, аб_тазик@ тай.ги 


Сегментация слабоконтрастных 
многопараметровых изображений 

по фазовым характеристикам 

метода ортогональной декомпозиции 


М.А. А1ехеует, Г.М. Уаонк, К.А. АКйтебйуп 
Манопа! Мттв Итуегуйу, Отртгоретоу5К, ИКгате 
49600 Риргорегоу5К, Кай Магх ау., 19, аб_тсавик @тай.ги 


Гом Соттая Мишрагатегет [тазе $евтетапоп оп Ва$е ой 
Рйазе Спагачетянсх ор Отйозопа! Ресотроятоп Метоа 


М.О. Алексеев, 1.М. Удовик, К.О. Ахметшин 
ДВНЗ «Нацщюнальний прничий ун1верситет», м. Дншропетровськ, Укра!на 
Укратна, 49005, м. Дншропетровськ, пр. К. Маркса, 19 


Сегментация низькоконтрастних багатопараметрових 
зображень за фазовими характеристиками методу 
ортогонально! декомпозици 


Рассмотрены информационные возможности нового метода анализа слабоконтрастных изображений. 
Суть метода основана на разделении анализируемых компонент многопараметрового (мультиспектрального) 
изображения на ортогональные составляющие, синтеза комплексных изображений и визуализации 
фазопространственных характеристик такого отображения. Представлены результаты реальной проверки 
работоспособности метода. 
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Введение 


Задача сегментации визуально неразличимых областей слабоконтрастных изоб- 
ражений (как обычных, так и многопараметровых) имеет большое прикладное значение 
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в таких областях, как медицинская диагностика, дистанционное зондирование земной 
поверхности и анализ изображений геофизических полей. В работе [1] были рассмо- 
трены информационные возможности метода ортогональной декомпозиции для реше- 
ния подобного класса задач, однако учитывались лишь амплитудные (яркостные) 
характеристики синтезированных изображений, а фазовая информация не учитывалась. 

Целью работы является демонстрация информационных возможностей фазовых 
характеристик метода ортогональной декомпозиции применительно к решению задач 
сегментации визуально неразличимых участков слабоконтрастных многопараметро- 
вых изображений. 


Ортогонализирующие преобразования 
многопараметровых изображений 


Для определенности первоначально рассмотрим случай, связанный с анализом 
многопараметрового (например, цветного) изображения { Г(х, у, п) }, где п — число 


параметров (диапазонов), соответствующих каждому пикселю с координатами (х, у). 
Если представить анализируемый ансамбль изображений {Т} в виде прямоугольной 
матрицы А размерностью (М * №)хи, где М и М - число пикселей вдоль осей хиу 


соответственно. Ортогонализация матрицы А возможна на основе использования раз- 
личных подходов [2]. В работе [1] нами были использованы два возможных подхода: 
метод сингулярного разложения и метод независимых компонент. 

В рамках первого подхода производится сингулярное разложение матрицы 24 
на основе выражения [3]: 


А=Ии\УТ, (1) 


где 0 = [#1 ОРвы О е — матрица сингулярных векторов размерностью ((М * №) хп), 
формирующих ортонормированный базис пространства, натянутого на столбцы мат- 
рицы А, и формирующая базис «собственных изображений» (СИ) [4] С(х, у, п); И - 
диагональная матрица, элементы которой и, (1=1,....п) являются сингулярными 
числами матрицы А. Они являются собственными значениями матрицы А'А, аУ- 
представляет собой ортогональную матрицу размерностью (пхп). 

В работе [1] нами анализировались СИ С(х, у, п) , которые представляют собой 
яркостные (т.е. амплитудные) характеристики синтезированного ансамбля С(х, у, п). 
Однако, учитывая ортонормированность С(х, у, п) , возможен синтез дополнительных 
трех фазовых характеристик для исходного цветного изображения: 


Ф, (х,у) = {[С, (х, у) + 16, (х,У)|, К =1,2,31 = 2,3;К =1. (2) 

Метод независимых компонент — это метод нахождения скрытых факторов или 
компонент в многомерных данных, особенность которого в том, что он позволяет 
выделять статистически независимые составляющие [5]. 

Основной принцип нахождения независимых компонент базируется на макси- 
мизации негауссовости анализируемого статистического распределения на основе 
использования отображения А => Р(х, у, п) ‚ где плотность распределения вероятностей 
Р4А(х,,х,...х, )} > РЕЕ)Р(Е, )...Р(Е,). В этом смысле метод независимых компонент 


обеспечивает статистическую, а не алгебраическую ортогонализацию. Но, тем не 
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менее, при его использовании так же возможен дополнительный синтез трех фазовых 
характеристик: 


Ч, (х,у)= Е. (х,у)+ Е. (х,У)|, К =1,2,31 = 2,3; К +1. (3) 


Однако, помимо метода сингулярного разложения, можно использовать еще 
один алгебраический метод ортогонализации на основе использования разложения 
Грама-Шмидта [6]. 

В качестве первого вектора выбирается у, =с,, где с,- первый столбец мат- 
рицы А, и затем каждый последующий вектор у,, ортогональный к векторам 


у, ,.... У, 1, рассчитывается по формуле: 


т 
_ _УнаС; 


о (4) 


... т 
и а 


Для любого номера { подпространство, порожденное исходными векторами 
с ›...с,, совпадает с подпространством, порожденным векторами уъ,....,у,. В резуль- 


тате векторы 4; =У, / |, ‚ которые можно развернуть в формате изображения, стано- 
вятся ортонормированными, что также позволяет синтезировать три фазопростран- 


ственные характеристики: 
О, (х,у) = {[0, (х, у) + 0, (х,У)|, К=1,2,31 = 2,3;К =1. (5) 
Таким образом, возможны три варианта синтеза фазопространственных харак- 
теристик при анализе слабоконтрастных многопараметровых изображений. 


Экспериментальные исследования 


Были исследованы три типа изображений: цветное космическое изображение; 
трехпараметровое медицинское изображение и «трехполевое» геофизическое изоб- 
ражение одного и того же участка местности. 

1. Космическое изображение. Конвертированное к градациям серого исходное 
цветное изображение представлено на рис. 1 а). Интерес в данном случае представляет 
анализ акватории морской поверхности с точки зрения обнаружения потенциальных 
залежей углеводородов. 

Исходное изображение (рис. 1 а) представляет собой классический пример слабо- 
контрастного изображения, поскольку перепады яркости на большинстве участков 
не превосходят 2%. Синтез нового изображения на основе трех СИ метода сингуляр- 
ного разложения (рис. 1 6) не решает задачи, и лишь использование фазо-простран- 
ственных характеристик Ф(х, у) (рис. 1 в) решает задачу сегментации. Сегментация 
изображения на основе фазовых характеристик Ф(х, у) метода независимых компо- 
нент (рис. 1 г) уступает в степени детализации методу сингулярного разложения, да 
и в программно-алгоритмическом аспекте этот метод представляется намного более 
трудоемким, по сравнению с методом сингулярного разложения. Отдельного внимания 
заслуживают фазопространственные характеристики О(х, у) метода ортогонализации 
Грамма-Шмидта (рис. 1 д), поскольку их сопоставление с фазовыми характеристиками 
метода сингулярного разложения (рис. 1 в) показывает, что они дополняют друг друга. 
Так, на рис. 1 д сегментируется нижний правый квадрат исходного изображения, а 
на рис. 1 в- три других квадрата. 
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Для сопоставимости результатов на рис. 1 е) представлен результат сегментации 
на 12 кластеров исходного изображения известным методом нечетких С-средних [7]. 
Сопоставление рис. 1 е) с рис. 1 в) и рис. 1 д) показывает, что этот метод проигрывает 
в степени детализации фазовым характеристикам метода ортогональной декомпози- 
ции, хотя если рассматривать ситуацию в целом, то результаты близки к друг другу, 
за одним исключением: затраты машинного времени на реализацию фазового метода 
ортогональной декомпозиции, для изображений размером порядка (300х 500) пикселей 
примерно на два порядка меньше, по сравнению с итерационным методом нечетких 
С-средних. 


Е д е 
Рисунок 1 — Космическое изображение: а — оригинал; 6 — синтез из трех собственных 
изображений метода сингулярного разложения; в — сегментация изображения на 
основе трех фазопространственных характеристик Ф(х, у); 

г_ сегментация изображения на основе трех фазопространственных характеристик 
Ч (х, у) метода независимых компонент; д — сегментация на основе трех 
фазопространственных характеристик ©(х, у) метода Грама-Шмидта; е — результат 
сегментации методом нечетких С-средних (12 кластеров) 

2. Медицинское радиологическое изображение. На рис. 2 представлено трех- 


параметровое медицинское ядерно-магнитно-резонансное томографическое изобра- 
жение опухоли головного мозга. 


{ 


а 6 В 
Рисунок 2 — Томограмма опухоли головного мозга: а — спин-решетчатая релаксация; 
6 — спин-спиновая релаксация; в — плотность протонов 
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Характерной особенностью данной ситуации является тот факт, что сама по себе 
опухоль хорошо видна, однако невидимой является область ее «скрытого влияния». 
Здесь сразу следует отметить, что ни использование метода нечетких С-средних, ни 
фазовых характеристик Ч(х,у) и О(х, у) не позволило решить поставленную задачу 
и лишь синтез фазопространственных характеристик Ф(х,у) метода сингулярного 
разложения позволил решить эту задачу (рис. 3). 


Рисунок 3 —Характеристики ортогональной декомпозиции (рис. 2) 
методом сингулярного разложения 


Из рассмотрения рис. 3 следует, что область скрытого влияния выделилась на 
фазопространственной характеристике Ф, (х, у) ‚ тогда как на яркостных характери- 
стиках С, область скрытого влияния не детектируется. 


3. Геофизическое изображение. На рис. 4 представлены изображения трех раз- 
личных физических полей одного и того же участка Земной поверхности. 


кь Ш 


Рисунок 4 — Изображения геофизических полей одного и того же участка 
Земной поверхности: а — гравитационное поле; 6 — магнитное поле; 
в — потенциальное электрическое поле 
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Основной задачей исследований геофизических полей является поиск аномаль- 
ных участков, с точки зрения обнаружения рудных залежей. При этом считается, что 
комплексирование различных геофизических методов позволяет повысить достовер- 
ность их обнаружения. 

На рис. 5 представлены результаты применения метода сингулярного разложе- 
ния для решения задачи сегментации таких потенциально интересных участков (все 
изображения цветные, но конвертированы к изображениям с градациями серого). 


а б Г 


Рисунок 5 — Конвертация к градациям серого цветных изображений: а — рис. 4; 
б-— С. › . метода сингулярного разложения; в — фазопространственные 
1,2,3 


характеристики Ф, › ‚; г- метод нечетких С-средних (12 кластеров) 


Из рассмотрения рис. 5 следует, что, как и в предыдущих случаях, использование 
фазопространственных характеристик метода сингулярного разложения позволило 
выделить несколько потенциально интересных участков, тогда как амплитудные 
(яркостные) характеристики сделать это не позволяют. 

Сопоставление рис. 5 в) с рис. 5 г) показывает, что степень детализации в фазо- 
вом методе ортогональной декомпозиции выше, по сравнению с методом нечетких 
С-средних, при том же выигрыше в машинном времени и программно-алгоритмической 
реализации. 


Заключение 


На основании проведенных исследований можно сделать следующие выводы. 

1. Для решения задач сегментации слабоконтрастных многопараметровых изо- 
бражений следует использовать фазопространственные характеристики метода син- 
гулярного разложения анализируемых ансамблей. 

2. Целесообразно сочетание ортогонализации слабоконтрастных многопараме- 
тровых изображений на основе комбинации метода сингулярного разложения и метода 
Грама-Шмидта. 

3. Сегментация слабоконтрастных изображений по фазовым характеристикам 
метода ортогональной декомпозиции примерно на два порядка является более быстро- 
действующей, по сравнению с методом сегментации на основе алгоритма нечетких 
С-средних, при несравненно более простой программной реализации. 
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ВКЕЗОМЕ 


М.А. Айехеуеу, 1.М. Удо, К.А. АКИтебруп 
Гом Соттая МиШрагатет [тазе Зеотетаноп оп Ва5е Рйа5е 


Сйагачетьунсу о7 Отповопай Ресотрозтоп Метоа 

То [1] сопяаегеа Фе шРюоптаНоп сара Ш@ез оЁ Фе тефо4 оЁ опвогопа| Чесотроз1оп 
(ю зоуе 1$ с1аз$ оЁ ргоетз, Баё юок шю ассоии оШу Фе атрШаае (Би тез$) 
свагацетзис$ оЁ Фе зупез17е4 птасе, ап рВазе шоптайоп — 1$ по (аКеп шю ассоипе. 

Тре ригрозе оЁ 15$ рарег 15 ю 4етопзиае Фе сБагасцет$ис$ оЁ рВазе шЮюппаНоп 
сара вез оЁ Ше плефо4 оЁ огфогопа! десотроз!оп у тезресе ю Ше {азК$ оЁ зестешайноп 
у1зиаПу шазипеизваЫе ю\/-сопёгаз( агеаз о пирагатаеег ппазез. 

То зоуе Фе ргоМепт$ оЁ зеотегианоп оЁР Фе 1о\-сопгаз птазез И пеед$ ю изе 
ши брагатеег рвазе-$расе сВагасепзИс$ оЁ фе тефоа оЁ зтеч]аг уаще десотрозвоп оЁ 
Фе апа[у7е4 ЪБапд$. 

Е 1$ аау1заЫе ю пих Фе огфогопайтайоп о? |ю\-сошгаз тшШарагатаег итахез Базе 
оп а сотЫпабоп оЁ зтечЧаг уае десотроз$!оп тео4 ап@ Фе Отат-ЗсВиа те@фод. 

Зестетаноп оЁ 1о\/-сопгазЕ птасез Бу Фе рБазе сБагацетзис$ ог Ше тефо4 оЁ 
огфогопа! десотрозвоп Бу абоие \о ог4ег$ оЁР таепииде 1$ тоге гезропяуе ш сотрайзоп 
У\и Фе тефо4 оЁ зеотещаноп Базе оп Фе С-теап$ Ёал77у аеопт а писб тпоге еазу 
зой\аге пиретешабоп. 


Статья поступила в редакцию 31.05.2012. 
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